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Abstract 

Penning trap measurements using mixed beams of 100 Mo - 100 Ru and Ge - Se 
have been utilized to determine the double-beta decay Q-values of 100 Mo and 76 Ge 
with uncertainties less than 200 eV. The value for 76 Ge, 2039.04(16) keV is in agree- 
ment with the published SMILETRAP value. The new value for 100 Mo, 3034.40(17) 
keV is 30 times more precise than the previous literature value, suffcient for the 
ongoing neutrinoless double-beta decay searches in 100 Mo. Moreover, the precise 
Q-value is used to calculate the phase-space integrals and the experimental nuclear 
matrix element of double-beta decay. 
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1 Introduction 

Neutrinos are one of the least understood fundamental particles. For half a 
century physicists thought that neutrinos, like photons, had no mass. But re- 
cent data from the neutrino oscillation experiments at SuperKamiokande pQ, 
SNO [2], and KamLAND [3j overturned this view and confirmed that the neu- 
trinos are massive particles. However, oscillation experiments can yield only 
the differences in the square of neutrino masses, therefore, no absolute mass 
values can be determined. In addition, another question remains concerning 
the fundamental character of neutrinos, whether being Dirac or Majorana 
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particles. Neutrinoless double-beta decay (O^/3/3-decay) is a process which can 
address both issues raised above. This decay process is forbidden according 
to the Standard Model of Particle Physics since it violates the lepton-number 
conservation and is only allowed if neutrinos are massive Majorana particles. 

The Heidelberg-Moscow Collaboration [4] has claimed the observation of the 
0^/3/?-decay using high-sensitivity [5.6J 76 Ge semiconductor detectors [7]. 100 Mo 
is another suitable nucleus to study the 0^/?/3-decay for the following reasons. 
100 Mo is rather easy to produce as enriched material and its 2z//3/3 matrix el- 
ement is known [8]. Furthermore, it has a high Q- value and the O^/5/3-decay 
rate scales with Q 5 . The MOON [9] and NEMO-3 [TO] experiments aim to use 
100 Mo for the O^/3/3-decay search. The signal for a O^/3/3-decay in a detector 
measuring the total energy of both electrons would be a peak at the position of 
the Q-value of the involved transition. Thus it is vital for the ongoing search 
of the 0^/9/?-decay that the Q- value is known to a fraction of the expected 
detector resolution. 

Additional motivation for the precise Q-value measurement is to calculate the 
precise phase-space integral G and the experimental nuclear matrix element of 
the two-neutrino double-deta decay [II] . The half-life Ti/ 2 of the double-beta 
decay process is expressed as 

[Tfa]- 1 = G 2v {M 2v ) 2 , (1) 



for the 2i//? / 9-decay and 

py a ]-i = G°»(M°») 2 ((m v )/m e ) 2 , (2) 



for the 0i//3/3-decay. Here M is the nuclear matrix element between the initial 
and final states of the decay, (m u ) is the effective neutrino mass [IT] and m e 
is the electron rest mass. For the 2i//9/9-decay case T 2 J 2 can be derived from 
experiment and G 2u calculated reliably. Hence using Eq. (1) the experimental 
M 2u is estimated. In the case of the Oz//?/?-decay the experimental half-life is 
unknown and hence we give the half-life of 100 Mo as a function of the effective 
neutrino mass for a given set of matrix elements using the G 0v value taken 
from this work. 

In this paper we present the precise double-beta decay Q-value of 100 Mo and 
report on the development of an off-line ion source for producing stable (or 
long-lived) mixed-ion beams for fundamental physics studies. 
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Fig. 1. (A): Schematic drawing of an off-line ion source to produce a variety of 
stable and long-lived mixed singly-charged ion beams. (B): Mass spectrum of ions 
ejected from the purification trap after buffer gas cooling for mass region A = 100. 
A gaussian fit to the data yields a FWHM of about 11 Hz, which corresponds to a 
mass resolving power of about 1 x 10 5 . 

2 Experimental method 



JYFLTRAP P2] is an ion trap experiment for high-precision mass measure- 
ments of radioactive ions ^3|14|15|ll6|17j . JYFLTRAP has been used for di- 
rect Qec~ value measurements of the superallowed beta emitters [T8fl9f20j and 
recently for double-beta decay. Our Q-va\ue measurement is based on a com- 
parison of the ion cyclotron frequency of a stored ion in a Penning trap |21j . 
The cyclotron frequency v c of an ion is given by 



1 9 R 
2ir m 



(3) 



where B is the magnetic field, m is the mass and q the charge of the ion. 
The ions of interest are produced at the Ion Guide Isotope Separator On-line 
(IGISOL) facility [22] which has an advantage that both mother and daughter 
nuclei can be produced in the same reaction. On the other hand, an off-line 
ion source was used to produce a mixed ion beam of long-lived (mother and 
daughter) nuclei. Thus it enabled the measurement of the cyclotron frequencies 
of the mother and daughter nucleus in a consecutive manner. Hence the Q- 
value can be obtained by using the following formula: 

Q = m m -m d = (— -l)m d , (4) 



where m m and are the masses of the mother and daughter ions and — is 
their cyclotron frequency ratio. The daughter nucleus was used as a reference 
ion and its mass excess value was obtained from ref. [23J. 
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In the case of the double-beta decay, 100 Mo (mother nucleus) and 100 Ru 
(daughter nucleus) ions were produced by using an off-line ion source shown 
in Fig. 1(A). The ion source consists of two electrodes and it is similar in size 
to the light ion guide used for indiced fusion reactions at the IGISOL facility 
[22]. One of the electrodes was designed to hold the metal plates or powder 
of the enriched element of need and the other one was the ground electrode. 
To create ions a continuous spark was ignited by applying 500 V between the 
two electrodes at a 15 mbar helium pressure. In general this off-line ion source 
can be used to produce any stable mixed-ion beam. In this work it was used 
at JYFLTRAP to create 76 Ge - 76 Se and 100 Mo - 100 Ru pairs of singly-charged 
ions. 

Ions were extracted from the off-line source by helium flow and guided by the 
sextupole ion guide (SPIG) into a differential pumping stage where they were 
accelerated to 30 keV and mass-separated with a 55° dipole magnet with a 
mass resolving power M/AM of ~ 500. The ions were then transported to a 
radiofrequency quadrupole (RFQ) structure where they were cooled, accumu- 
lated and bunched [24J. Finally, they were injected into the double Penning 
trap system. 

The JYFLTRAP Penning traps are placed inside the warm bore of a 7 T 
superconducting magnet and they are separated by a 2-mm channel. The first 
trap is called the purification Penning trap, where the mass-selective buffer- 
gas cooling technique is applied for further axial cooling and isobaric cleaning 
[25] . The mass resolving power of the purification trap was on the order of 10 5 . 
Figure 1(B) displays a mass spectrum of ions ejected from the purification 
trap after the buffer gas cooling. The purification trap was operated to center 
selectively the mass region A = 100. The one of the selected ion samples was 
transported to the second trap called the precision trap. 

In addition, time-separated (Ramsey method) [26] dipolar excitation was ap- 
plied in the precision trap to ensure a single ion species is selected [27]. This 
was performed in the following way: at first the mass-selective reduced cy- 
clotron frequency u + (removal frequency) was applied for one ion species as 
two time-separated fringes of 5 ms with a waiting time of 20 ms (5-20-5 ms). 
This increases the reduced cyclotron radius of the unwanted ions. At the end 
of the excitation the ions were sent back to the purification trap, thus the 
excited ions (i.e. the unwanted ions) will not pass the 2-mm channel between 
the traps. This method allows us to clean even isomeric states that are only 
200 keV higher than the ground state (in the case of 54 Co) [28]. In the purifi- 
cation trap, the buffer-gas cooling technique was re-applied and finally a pure 
ion sample was transported to the precision trap for the cyclotron frequency 
(u c ) measurement. 
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At first, the cyclotron frequency was determined by employing the time-of- 
flight (TOF) technique [29] with an excitation time of 100 ms in a single- 
fringe quadrupolar excitation scheme (conventional excitation scheme). Once 
the cyclotron frequency was determined, a quadrupolar Ramsey excitation was 
applied in order to determine the cyclotron frequency with higher precision 



A typical TOF resonance is shown in Fig. 2 (top panel) for 100 Ru ions with a 
conventional excitation scheme and an excitation time of 100 ms. The middle 
and bottom panels in Fig. 2 show the Ramsey resonances for 100 Ru and 100 Mo 
ions, respectively. A Ramsey excitation with two time-separated fringes of 
50 ms and a waiting time of 400 ms (50-400-50 ms) was applied for these 
particular cases. 



3 Analysis and Results 

The statistical standard deviation cr(v c ) of the cyclotron frequency u c can be 
estimated as [3T] 



where iV = number of detected ions, T ex = excitation time and K is an empiri- 
cal constant independent of ion number and excitation time. In this experiment 
the precision was maximized by collecting a large number of ions and by ap- 
plying long excitation times. Also a higher v c improves the relative precision 
achieved with a strong magnetic field (see Eq. (1)). In addition, the Ramsey 
excitation yields a narrower central peak compared to the conventional excita- 
tion scheme. Therefore, using the Ramsey excitation the precision is improved 
by a factor of about 2 or more [32]. All these factors together allow us to 
determine the Q-value of the double-beta decay of 100 Mo with high-precision. 

For the consistency of the setup we have measured the Q-value of the double- 
beta decay of 76 Ge which was measured in 2001 at SMILETRAP [33J and 
recently again in 2007 applying the Ramsey excitation method [31]. Our value 
agrees with the SMILETRAP value (2039.006(50) keV) [33] within the error 
bars. For 76 Ge 25 resonances were collected with 76 Se as a reference. A Ramsey 
excitation scheme was applied with two time-separated fringes of 25 ms with 
a waiting time of 250 ms (25-250-25 ms). For 100 Mo a total of 68 resonances 
were collected with 100 Ru reference. These were grouped in two sets with 
Ramsey excitation times of 25-250-25 ms and 50-400-50 ms. Figure 3 shows 
the scattering of the individual frequency ratios r« relative to the weighted 
average value of the frequency ratio f = 1.000 032 6069(18). The inner and 
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Fig. 2. Time-of-flight (TOF) resonances of 100 Ru and 100 Mo ions from the precision 
trap. The top panel represents a typical TOF resonance using the conventional 
excitation scheme with an excitation time of 100 ms. The middle and bottom panels 
display the TOF resonances for 100 Ru and 100 Mo ions using the Ramsey excitation 
pattern with 50-400-50 ms, respectively. 



outer statistical uncertainties (Sf) [35] of the weighted average frequency ratio 
are 1.6xl0 -9 and 1.8 xlO -9 , respectively. The ratio of these values (Birge 
ratio) is very close to 1 which confirms that the scattering of the data is 
statistical. However, the larger error (in this case outer error) was used to 
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Fig. 3. Measured cyclotron frequency ratios fj between Mo + and Ru + relative 
to the weighted average value f. 

determine the final uncertainty. 

In the determination of the cyclotron frequency ratio the following systematic 
uncertainties were taken into account: The number of ions present in the 
trap can cause a shift in the cyclotron frequency. This is taken into account 
by plotting the center value of the cyclotron resonance as a function of the 
detected number of ions [36]. The cyclotron frequency equivalent to one ion 
in the trap was determined from this plot via a linear extrapolation to the 
observed 0.6 ions (detector efficiency = 60%). This extrapolated cyclotron 
frequency was taken as a final value. Hence the uncertainty due to the count 
rate comes together with the statistical uncertainty. 

The drift of the magnetic field was taken into account by an interpolation of 
the reference frequencies measured before and after the cyclotron frequency 
measurement of the ion of interest. This linear interpolation does not take into 
account the short term fluctuations of the magnetic field. This was taken into 
account by adding AB/B = 3.22(16)x 10 _11 /min [TB] multiplied by the time 
difference between the two consecutive reference measurements quadratically 
to the uncertainty of the frequency ratio. As the Q- value is determined by the 
cyclotron frequency ratio between the mother and daughter having same ^, 
mass-dependent and other systematic uncertainties (if it exists) should cancel. 

The Q-value can be derived from the final weighted average frequency ratio 
using Eq. (4). The term inside the parenthesis in Eq. (4) is small (~ 2 — 3 x 
10~ 5 ), thus the uncertainty contribution from ma to the Q- value is negligible. 
Figure 4 shows the weighted average Q- values of 100 Mo with their uncertainties 
for two different excitation times and these two sets of data agree each other. 
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Table 1 

Weighted average frequency ratios (f) and the Q-values of 76 Ge and 100 Mo mea- 
sured at JYFLTRAP. The final uncertainty normalized with the x 2 1S given in the 
parenthesis. # and T ex represent the number of doublet measurements and exci- 
tation scheme, respectively. Final frequency ratios and Q-values are indicated in 
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Fig. 4. Q-values obtained from individual measurements of Mo in June 2006 by 
using 100 Ru as a reference. Triangles and circles correspond to different excitation 
times (see legend). The data collection was stopped every night before the automatic 
Bq dump at 3:00 AM. The statistical, count-rate and magnetic-field fluctuation 
uncertainties have been taken into account for each set. 

The final weighted average cyclotron frequency ratios f and Q-values of 76 Ge 
and 100 Mo with their corresponding uncertainties in the parenthesis are given 
in Table 1. 



4 Discussion 



The precise 100 Mo Q-value is used to calculate the phase-space integrals G of 
the 2i//3/3-decay and 0z//3/3-decays. A detailed formula for these calculations can 



8 



Table 2 

Phase-space integrals for 100 Mo. The axial- vector coupling constant g^ = 1.0 and 
1.254 and electron rest mass m e = 510.998918 keV were used in the calculations. 
The uncertainty for G is estimated solely from the Q-value uncertainty. 



Q-value / keV 


Q2u x 1Q -18 


for 


G 0u x 10 -u 


for 




9A = 1.0 


9A = 1-254 


9A = 1.0 


9A = 1-254 


3034.40(17) 


3.8179(19) 


9.4409(47) 


1.8915(4) 


4.6772(10) 


3035(6) 


3.826(67) 


9.46(17) 


1.893(14) 


4.681(34) 



Table 3 

A comparison between the calculated and experimental (Ex) nuclear matrix ele- 
ments M 2iy for 100 Mo. For detailed notation see ref. [llj . 



Ex 


Ex 


QPRA 


QPRA 


SPRA 


SU(3) 


SU(3) 


for 


for 


EMP 


EMP 


WS 


SPH 


DEF 


<M=l-0 


<M=1.254 


m 


m 


@QJ 


m 


m 


0.192(4) 


0.122(4) 


0.256 


0.197 


0.059 


0.152 


0.108 



be found in ref. [TT]. Results of the calculations are summarized in Table 2. The 
uncertainty in G is estimated solely from our Q-value uncertainty. The other 
uncertainty (in addition to the uncertainty from the Q-value) in G comes from 
the Fermi function approximation for the exact solution of the Dirac equation 
for a homogeneously charged sphere. The relativistic Fermi approximation 
limits the accuracy of G to some three digits. 

Using Eq. (1) the value of the nuclear matrix element M 2v is estimated exper- 
imentally to be 0.122±0.004 where the precise phase-space integral value is 
taken from this work and the recommended half-life Tu2 = (7.1 ± 0.4) x 10 18 
years is taken from [37J. A comparison between the computed nuclear matrix 
elements and the experimental one is shown in Table 3. Theoretical values of 
M 2u vary by a factor of five whereas the experimental value lies in between. 
The matrix elements 0.197 and 0.152 of the spherical QRPA and SU(3) theo- 
ries are consistent with the range of the experimental matrix elements. On the 
other hand the matrix element 0.108 computed by using the deformed SU(3) 
theory falls out side this range. This would suggest that deformation is not 
very important for the ground-state to ground-state 2v(3f3 decay of 100 Mo. 

In the case of 0i//3/3-decay finding an estimate for the experimental value of the 
nuclear matrix element M 0u is not so straight-forward since the decay half- 
life is unknown in Eq. (2). However the expected half-life can be plotted as a 
function of the effective neutrino mass as shown in Fig. 5. We have used our 
value for the phase-space integral and a set of nuclear matrix elements M 0u 
calculated by using the pnQRPA }42ll43j and RQRPA [13] models. In these 
models the nucleon-nucleon short-range correlations have been accounted for 
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Fig. 5. Expected half-life band as a function of the effective neutrino mass for Mo. 
The ranges of the nuclear matrix element M° u come from |42j in the upper panel 
and from [43j in the lower panel. The values of the phase-space integrals were used 
from this work. 

by two different methods, the Unitary Correlation Operator Method (UCOM) 
and Jastrow. For both models, the pnQRPA and the RQRPA, the range of 
computed matrix elements stems from the uncertainty in the value of the 
axial-vector coupling coefficient, = 1.0 - 1.25, and from the variation in the 
value of the proton-neutron particle-particle interaction strength g pp , used to 
fit the experimental range of the 2vf3f3 half-life. 
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5 Conclusion 



In this article an off-line ion source is demonstrated for creating stable (or 
long-lived) mixed beams of single- charged ions. A precise double-beta decay 
Q-value, 3034.40±0.17 keV for 100 Mo, is determined, a candidate for searching 
the evidence of 0z//3/3-decay. The Q-value precision of about 200 eV was ob- 
tained for mass 100 region. Experimental phase-space integral G is obtained 
and using this value an experimental M 2v value is derived for 100 Mo. In the 
case of neutrinoless double-beta decay half-life is estimated as a function of 
the effective neutrino mass for a set of nuclear matrix elements. 
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